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TECHNICAL NOTE 2321 

ANALYSIS OF TEMPERATURE DISTRIBUTION IN LIQUID -COOLED TURBINE BLADES 
By John N. B. Livingood and W. Byron Brown 


SUMMARY 

The temperature distribution in liquid-cooled turbine blades 
determines the amount of cooling required to reduce the blade tempera- 
ture to permissible values at specified locations. This report pre- 
sents analytical methods for computing temperature distributions in 
liquid-cooled turbine blades, or in simplified shapes used to approxi- 
mate sections of the blades. 

The individual analyses are first presented in terms of their 
mathematical development. By means of numerical examples, comparisons 
are made between simplified and more complete solutions and the effects 
of several variables are examined. Nondimens ional charts to simplify 
some temperature-distribution calculations are also given. 

It was found that for blade materials having low thermal conduc- 
tivities (corresponding to blade materials in current use) a one- 
dimensional spanwise temperature -distribution equation is applicable 
near a coolant passage. The cooled part of the blade, irrespective of 
blade length, was found to remain at a determinable uniform temperature 
called the prevalent blade temperature. The prevalent blade tempera- 
ture increased about 150° F for an increase in the difference between 
the effective gas temperature and the coolant temperature of 1000° F 
for a range of effective gas temperature from 2000° to 5000° F. 

In cases where rim cooling is insignificant, it was found that a 
one -dimens ional chordwise temperature -distribution equation was suf- 
ficiently near the relaxation solution for the actual blade shape for 
a first-order approximation. 

INTRODUCTION 

I 

A limitation on design and performance of aircraft gas-turbine 
power plants is the. strength of the turbine materials, which decreases 
as temperature increases. This limitation may be greatly alleviated, 
even when nonstrategic materials are used, by the application of tur- 
bine cooling. An extensive study of both liquid and" air cooling of 
turbine blades has been carried out at the MCA Lewis laboratory 
since 1945. 
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The temperature distribution in a turbine blade determines the 
amount of cooling required to reduce the blade temperature to permis- 
sible values at specified locations. Analytical studies have been con- 
ducted to predict temperature distributions in cooled turbine blades. 

A summary of the analytical methods developed for calculating spanwise 
temperature distributions for three types of air-cooled turbine blade is 
presented in reference 1. The present report summarizes the analytical 
methods developed for calculating temperature distributions in liquid- ' 
cooled turbine blades. Parts of these analyses are presented in refer- 
ences 2 and 3. 

Because the calculation of a generalized three-dimensional temper- 
ature distribution through a liquid-cooled turbine blade is too tedious 
to be of practical value and because some knowledge of the distribution 
in a radial direction is necessary, the first investigation presented 
herein is limited to a one -dimensional spanwise distribution. This 
distribution, however, is only valid in the immediate vicinity of the 
coolant passages. In order to study more accurately the temperature 
distribution in the trailing section of the blade, a three-dimensional 
solution is determined for a rectangular parallelepiped used to approx- 
imate a blade trailing section. This solution serves to determine the 
importance of blade thermal conductivity and of distance from the coolant 
passage. Because numerical calculations for a rectangular parallele- 
piped produce a constant spanwise temperature in the region of the blade 
beyond the influence of rim cooling (that is, in the region of the blade 
unaffected by conduction to the rim), one -dimensional chordwise tempera- 
ture distributions are then calculated for shapes approximating various 
parts of a blade cross section) these investigations are valid for the 
part of the blade beyond the influence of rim cooling. Finally, a two- 
dimensional temperature distribution through the cross section of an 
actual blade shape is obtained in order to determine the accuracy of 
the previously deter min ed approximate solutions and to study the effect 
on t em perature distribution of a peripheral variation of gas -to-b lade 
heat -transfer coefficient . 

Numerical examples based on coefficients available from unclassi- 
fied sources are included. 


METHODS OF ANALYSIS 

One -dimensional spanwise, three-dimensional, one-dimensional 
chordwise, and two-dimensional chordwise temperature -distribution equa- 
tions for liquid-cooled turbine blades or for simplified shapes used to 
approximate sections ©f liquid-cooled turbine blades are presented. 
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For simplicity in presentation, the individual analyses are first formu- 
lated in terms of the mathematical methods. The comparative applica- 
bilities of the methods of analysis are established in later sections 
of the report. . 

For all the analyses presented, the following conditions are 
assumed: 

( 1) The coolant temperature is constant at the average temperature 
of the coolant, or the coolant forms a constant -temperature reservoir. 
This assumption is valid when the change in coolant temperature is small 
in compaxison with the temperature difference between the gas and the 
coolant. 

(2) The blade -to-coolant heat-transfer coefficient is constant. 

(3) The thermal conductivity of the blade is constant. 

, - (4) The effective gas temperature at all parts of the blade is the 
same. 


(5) Radiation effects axe considered to be included in the heat- 
transfer coefficients. The following individual analyses axe presented: 

One -dimensional spanwise temperature distribution. - A one- 
dimensional spanwise calculation gives a radial temperature distribution 
valid in a part of the blade near a coolant passage. In this particular 
case, the radial temperature distribution is carried through the blade 
and through the turbine rotor. Heat transfer between the turbine rotor- 
and the fluid on the outside of the rotor is assumed to take place and 
a constant value of gas-to-blade heat-transfer coefficient is assumed. 
For simplicity, a blade of constant cross-sectional area and perimeter 
is considered. 

For blades with little taper, the method is applicable if average 
cross-sectional area and perimeter are considered. In general, however, 
a numerical solution is necessary for tapered blades. 

Three-dimensional temperature distribution through simplified 
trailing section. - A three-dimensional temperature distribution for the 
trailing section - of the blade is obtained by considering this part of 
the blade as a rectangular parallelepiped. The blade -root temperature 
is considered constant, the gas-to-blade heat-transfer coefficient is 
constant, and it is assumed that no heat flows across the median plane 
of the simplified blade section. 
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An approximate solution for tapered blades is obtainable by -this 
method if the blade is considered in small sections, average values of 
cross-sectional area and perimeter are considered for each section, and 
appropriate compatability relations are satisfied. 

One -dimensional chordvise temperature distributions through simpli- 
fied shapes. - Rectangular and trapezoidal sections, concentric circle 
annuli, and sections between parallel plates are considered as the sim- 
plified shapes vised to approximate parts of a turbine -blade cross sec- 
tion. The gas-to-blade heat -transfer coefficient is considered constant. 
For comparison purposes, the rectangular and trapezoidal sections are 
constructed so that lengths and areas are identical. 

Two-dimensional chordvise temperature distribution. - Numerical 
two-dimensional temperature distributions through the cross section of 
an actual blade shape are deter min ed - Gonstant and variable gas-to- 
blade heat-transfer coefficients are considered. 


One -Dimensional Spanwise Temperature Distribution 

The spanwise temperature -distribution equations for a liquid- 
cooled turbine, a section of which is shown in figure 1, are derived in ' 
reference 2 and are reproduced in appendix B. In this spanwise case, 
the investigation was extended from blade tip to rotor hub . 

For convenience, the turbine was divided into four sections 
(fig. 1), and the temperature-distribution equations were obtained from 
heat balances for differential elements in each section. Inasmuch as 
chordvise conduction was neglected, the following equations, valid only 
in the neighborhood of the coolant passages, were obtained. For the 
uncooled section of the blade (fig. 1, section 1), it was found that 

Tg,i = ^g,e ~ ^1 cos k -s^( x i + Cg) (1) 

where C-j_ and Cg are integration constants and 


1 



(All symbols are defined in appendix A.) The value of C2 is deter- 
mined by use of the boundary condition at the blade tip (where heat 
enters by convection only) . 
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/ar B \ 

" k B A B,l \d^j T = Mb, 1 ( T g,e- T B,T) 


or 


tanh 


= ^o_ = ( Mb,! 
\Mb > 


1 

2 


For the cooled section of the "blade (fig. 1, section 2), it was 
found that 


T B)2 = T B,p “ C 3 e “ C 4 e 


0Sx 2 _ -#x 2 

A 

where C3 and C4 are integration constants. 


' (2) 



Mo T g,e + h i,2 7 i,2 T Z 

Mo + h i,2 Z i,2 


( 2 a) 


and 

1 

(Mo + h i,2 2 i,2\ 2 

*- X %a B;2 -) 

The prevalent blade temperature T-g^p is the temperature the blade 
would assume if no heat flows* in or out of the blade ends. 

For the rim section (fig. 1, section 3), with assumed constant 

.area. 


L rim 


$ X3 -<?X3 

= $4- Cce + C^e 


(3) 


where C 5 and C 0 are integration constants. 
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^ Jir 5jav* I a + ( h i,5 Z i,3)av 


^rim Arirn 


and T a is the fluid temperature on the outside of the rotor and 
( h i ; 3 l i,3^av is the avera S e of the values of ^ 3 found for the 

blade coolant passages, the rim circumferential passages, and the radial 
passages through the rotor (fig. 1) . 


For the rotor section (fig. 1, section 4), with assumed constant 
strength, 

T r = /+ C ? cosh^. 

where Gj is an integration constant, and where 


(4) 


£ _ 4=nr 4,av ll a T a + (kj,4^i,4)av T 2 
^ rtr 4,av^ L a + ^i,4^i,4^av 


and 




4nr 


4,av^a + ^i,4^i,4^ av 


^ A R 


v 


Values of the six integration constants C^, C 3 , C 4 , C5, Cg, 

and Cy are found by solving simultaneously the six equations resulting 
from equating temperatures and heat flows at the various junction points 
of the four sections of the turbine. (Subscripts a and p used after 
numerical subscripts designate the end nearest the blade tip and the end 
nearest the rotor hub of the sections to which the numerical subscripts 
refer, respectively.) At the junctions of sections 1 and 2, 


T g,e - C 1 coshjrf(x 1#p + C 2 ) = T B , p - C 3 e - C 4 e"^ 2,a 


(5) 


k B A B,l,P C l‘ s/sinll ^ x l,P + C 2^ = ^ k B A B,2,a 


#x -<2x . 

C,e 2 > a - C A e 2 ' a ) + 


/ #x -^x 9 \ 

h i,l( A B,l,P - A B,2,a) K,p - C 3 e 2 ' a - C 4 e - T z) ' < 6 ) 
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Equation (6) equates the heat leaving section 1 to the sum of the heat 
entering the metal of section 2 and the heat entering the coolant at 
that part of the blade where the inlet and outlet passages are con- 
nected. An approximation in the procedure has been introduced at this 
point because a separate blade section for the part of the blade con- 
taining the connecting passage between inlet and outlet passages is not 
introduced. A numerical calculation showed that use of such a section 
would slightly decrease the temperature at the blade tip. 


At the junction of sections 2 and 3, 


T B,p " C 3 e 




<?x, 


- c 4 e 


C 5 e 


+ C 6 e 


3, a 


( 7 ) • 


L g,e-%p + C 3 e 


, #x 2p -*X2,p\ 

VB,2 # ( C 3 e C 4 e ~ ) 

^,2 hA ' 2 V 

A r im / *” ry *** fvN 

“ **B*B,Z (^6 e 3 ’ a - C 5 e 3 ’ a ) 




-<#x, 


+ C 4 e 


2,P 


( 8 ) 


Equation (8) equates the sum of the heat leaving section 2 and that 
entering section 3 directly from the hot gases to the total heat entering 
section 3. 


At the junctions of sections 3 and 4, 


$x, 


3,P 


- <£x. 


+ Cj-e + C g e 3 ' P = f + C ? cosh 


. O) 


£x. 


k B *(- C R e 


: 3,P - ^ X 3,P\ A : 


+ c 6 e 


\ A-o 

J - ^ 7 — %jfc 7 Sinhjfr a (10) 


Three-Dimensional Temperature Distribution through 
Simplified Trailing Section 

The spanwise temperature-distribution equations previously pre- 
sented are valid only in the neighborhood of the coolant passages-. Most 
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currently used turbine blades are so shaped that coolant passages can- 
not he located near the trailing edge. Because of the .impossibility of 
placing reasonably large coolant passages near the trailing edges of 
conventional turbine blades and because the trailing sections seem most 
likely to be the hottest parts of the blade, other detailed studies were 
made of temperature distributions . The first study was devoted to 
approximating the trailing section by a rectangular parallelepiped and 
to determining the three-dimensional temperature distribution through 
this parallelepiped (fig. 2). 

The median plane of the rectangular parallelepiped was chosen as 
the plane z = 0, and symmetry was assumed about this plane (no heat 
flow was assumed to cross this plane, as previously stated). Boundary 
conditions at the blade trailing edge and blade tip were simplified by 
assuming the blade width and blade length to be extended by a distance 
equal to one -half the blade thickness. In figure 3, the edge MN at 
temperature T gained some heat from the gas stream. The extended 
surfaces NO and MP were at nearly the same temperature T and no 
heat entered the edge OP; these surfaces therefore gained practically 
the same amount of heat as the actual exposed edge. The validity of 
this approximation is discussed in reference 4 (pp. 216 and 217). Dis- 
tances increased by t/ 2 are denoted by primes. 

The derivation of the three-dimensional temperature distribution, 
originally derived in reference 3, is reproduced in appendix. C. The 
differential equation, in final form, for the temperature distribution 
is found to be 


d 2 e d 2 e 

dx' 2 V 2 dz 2 


(ID 


where x' and y' denote x and y increased by t/2, respectively. 
A solution of equation (ll) satisfying the boundary condition’s at 
x 1 s y 1 s z = 0 is 

00 00 

(•#m,n cos y n x' coshi y 1 cos^z + 

D3f= 1 n=l 7 

0 m,n cosh ^m,n x ' cos cos ^m z ) ( 12 ) 

where jC y i£, Jl, Jf, 0, 9P, and 91 are constants. Relations among 
the constants are 
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and 


tan 


/ 1 = 
m 2 


k B ^ r 


m 


h T 

o 

2k B 

T 

“ 


(18) 


It is at once obvious that from equation (16) many values of 
result. Because of the periodicity of the tangent function, equa- 
tions (17) and (18) likewise have many solutions. 

Values for all the constants in equation (12) can thus be found 
and the temperature can now be computed at any point in the rectangular 
parallelepiped. In reference 3, appendix B, it is shown that suffi- 
ciently accurate results can be obtained by using m = 1 only. 


One -Dimensional Chordwise Temperature Distributions 
through Simplified Shapes 

Because the spanwise temperature distributions are valid only in 
the neighborhood of the coolant passages and because the three- 
dimensional approximate solution resulted in a constant spanwise blade 
temperature in the part of the blade beyond the influence Of rim cool- 
ing, one -dimensional chordwise temperature distributions were deter- 
mined for sections of a liquid-cooled turbine blade that can be approx- 
imated by simple shapes (fig. 4) . Rectangular and trapezoidal approx- 
imations for blade trailing sections were considered. Trapezoidal 
sections may also be used to approximate the leading section in some 
blades. In addition, analyses were made for leading sections approx- 
imated by concentric circle annuli and for the sections of blades with 
very little metal between blade outer surface and coolant passages 
approximated as regions between parallel plates. The temperature dis- 
tributions through these simplified shapes, derived in appendix D and 
valid in regions beyond the influence of rim cooling, follow. 

Rectangular trailing section of blade removed from influence of 
rim cooling. - For a blade trailing section approximated by a rectangle 
(fig. 4(a)) the temperature-distribution equation is 

h * 

kj (T g^ e " T ^ cosh 

6 — ■ - -■■■ ■■ , ■■ — ■■■ ■ — 

h i 

cp sinh cp j ' + t — cosh cp j ’ 
k B 


( 19 ) 
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where 


Q ~ T g,e-TB 


and 


cp = 


.V 


1 

2 


Trapezoidal trailing section of blade removed from influence of 
rim cooling* - For a trapezoidal approximation to the blade trailing 
section (fig. 4(b) ), the temperature distribution is 


( T g,e- T I> [%(i£:L>J 0 <i5> + iJiUyiHodti] 


e = 




[i^dyy iy] - [yiyi^dy] ♦ ^ r 


( 20 ) 


where 


£ = 2K 


y* + 


T.U - tan +) 


2 tan i|r J 
1 


5 = 2K 


Tl (l - tan i|r) 


L 2 tan \|r J 


? 2 = 2 k 


j' + 


T 1 (l - tan \|r) 


2 tan i | r J 

1 

v2 


K = 


\k B sin i| i / 


T -T 

, -1 2 1 
f = tan . — K-r- 

2j 


r = %(!£•,_) J 0 (i£ 2 ) + UxCiyiHoCi^) 
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and J Q , iJj^, H Q , and iHj are Bessel functions. For a wedge- 
shaped .section the temperature -distribution equation (20) reduces to 



2K 2 ko 

W 

2K 2 kB 




( 21 ) 


Concentric -circle -annulus approximation for blade leading section. - 
The equations expressing the blade -metal temperature at the coolant 
passage wall Tg^ ^ > at the leading edge of the blade Tg^ 0 , and at any 

point in the annulus Tg are found to be (fig. 4(c)) 


and 


T B,i = + 


Ttj ^ = 


2k B 


(V 0 \ 

D i h i log e\p^y 


( T B,o _T B,i) 


vr 


J 



) ( T B,o' T B,i) 


(22) 


(23) 


log- 


T B “ T B,i + 


/ _D \ 

v>J 


l0g e Id7 


D, 


( T B,o- T B,i) 


(24) 


where 


^B,o _T B,i 


1 + 



(25) 


Section of blade approximated by parallel plates. - The equations 
expressing the blade-metal temperatures at the coolant-passage wall 
(fig. 4(d)) Tg ± and at the blade outer edge Tg^ 0 are 
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T : + (snr) (h ° T g,e+*iTj) 

m ' O i/ 

B>1 ” / k* \ 

(26) 


1 ♦ Ua) <h ° +h i> 


and 

T g>e + (qh h.) ^ o T g,e + iV 
m \ O 1/ 

2,0 ~ / kn \ 

1 + (ohjhj <VA> 

(27) 


where o is the distance between the parallel plates. 


Two-Dimensional Chordwise Temperature Distribution 


The temperature at any point on the blade cross section may be 
found by solving the Laplace differential equation given in terms of the 
temperature difference 6 


9 a 2 9 

dx 2 5 y 2 


0 


(28) 


where X and Y are the. Cartesian coordinates in the blade cross sec- 
tion. The boundary conditions, expressed in terms of the partial deriv- 
ative of 6 in the direction normal to the. boundary are 


3e h o 
35“ % 6 


at the outer boundary and 


30 

5n 


h. 


= kg C^e-V 6 ) 


at the coolant -passage boundary. 


(29) 


(30) 


A closed solution to equation (28)' cannot be obtained because of 
the impossibility of applying the given boundary conditions along the 
odd-shaped boundary of a turbine blade (fig* 5). A numerical solution 
is available, however, by application of the relaxation method (refer- 
ence 5) • A sketch of the blade cross section is covered by a network 
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of points (a square network was chosen) . Large net spacings are recom- 
mended at first, and the insertion of additional net points may be made 
at any time during the solution, thus permitting the use of a final 
network of any desired size. Temperature estimates for the net points 
can be obtained from the solutions of the chordwise blade -temperature - 
distribution equations for approximate shapes previously presented and 
from application of the boundary conditions given by equations (29) and 
(30). Residuals, which may be considered as interior heat sinks, can 
then be calculated at each net point whose immediately neighboring 
points remain within the boundary from the relation 

Qq .= e-L+02+03+04-400 (31) 

where the subscript zero denotes a point in the blade cross section and 
subscripts 1, 2, 3, and 4 denote the points in the square surrounding 
the point with zero subscript. 

The object of the relaxation is to reduce the values of Q to 
zero, or as close to zero as possible. When the relaxation equa- 
tion (31) is_employed,_the following procedure is used. A change in 
Qq alters Q-[_, QgJ %■> ^4 ^y th® same change and Qq by minus 

four times this change, all other values of 9 remaining fixed. 

Equation (3l) is the finite -difference equation corresponding to 
the partial differential equation (28). 'For net points some of whose 
immediately neighboring points lie outside the boundary, equation (3l) 
must be modified. For example, if point 1 lies outside the boundary 
and point 5 is the boundary point between 0 and 1, the following equa- 
tion applies : 

% - w e 4 + t - ( 3 + IH < 32 > 

where d is the ratio of the distance between points 0 and 5 to the 
net spacing. 

Corresponding changes in the values of Q result from the use of 
equation (32). No harm is done by overestimating the final values of 
6 y as successive calculations will establish them again. Continued 
relaxation eventually reduces all the residuals as desired and the blade 
temperatures can finally be obtained from the definition of Q, 


0 = T g,e-*B 
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APPLICATION OP ANALYSIS 

Typical numerical examples for a sample blade are presented to 
illustrate temperature trends and to determine the effects of various 
factors on the temperature distributions in a liquid-cooled turbine 
blade. A turbine blade whose external shape is similar to that of the 
rotor -root section of a conventional gas-turbine design was selected 
as the sample blade. Two 0.25-inch-diameter coolant passages were 
assumed in the blade (fig. 5), connected near the blade tip by a cross- 
over passage, and the cooling was assumed to occur through forced con- 
vection. No allowance was made in these calculations for the effects 
of free -convection currents that might be present. A gas. flow of 
55 pounds per second (equivalent' to a mass velocity of about 58 lb/ 
(sec)(sq ft)) and a coolant flow of 7 pounds per minute per blade 
were assumed (to insure turbulent flow and to remove the dependence of 
the blade-to-coolant heat-transfer coefficient on the length-to- 
diameter ratio of the coolant passage) . 

For the one -dimensional spanwise case, in which the cooling was 
carried to the rotor hub, two 0.25 -inch-diameter coolant passages 
running circumferentially through the rim and ten 0.50-inch-diameter 
coolant passages running radially through the rotor were assumed 
(fig. l) . For this case, water, ethylene glycol, and kerosene were all 
considered as possible coolants. Other calculations were made only for 
water as the coolant because of its superiority over the other coolants 
considered. 

Blade-to-coolant average heat-transfer coefficients were calculated 
by use of formulas presented in reference 6, page 168 for turbulent flow 
and page 190 for laminar flow (for ethylene glycol only). Gas-to-blade 
average heat -transfer coefficients were calculated by use of formulas 
given on page 236 of reference 6. The following coefficients were used 
in the analysis: 


Outside heat -transfer coefficient, h Q 

Btu/(hr) (sq ft) (°F) 222 

Btu/(sec) (sq ft) (°F) 0.06167 

Inside heat -transfer coefficient, h^ g 

Water, Btu/(hr) (sq ft) (°F) 2370 

Btu/(sec) (sq ft) (° f) 0.6583 

Ethylene glycol, Btu/(hr) (sq ft) (°F) 649 

Btu/(sec)(sq ft)(°F) . . . . 0.1803 

Kerosene, Btu/(hr) (sq ft) (°F) 510 

Btu/(sec) (sq ft) (°F) 0.1417 
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Other numerical values used in the calculations were (figs. 1 and 
4): 


Effective gas temperature, Tg^ e .> °F 2000 - 5000 

Average coolant temperature, T^, °F . . 200 

Number of blades, Z 55 

Area of blade cross section 1, Ag, sq ft 0.00198 

Area of blade cross section 2, Ag, ' sq ft . . '. 0.0013 

Area of cross section 3, A3, sq ft . . 0.312 

Area of cross section 4, A4., sq ft 0.236 

Blade outside perimeter, Z 0 , ft 0.2542 

Blade inside perimeter, ^i,2-> ^t 0.131 

Average radial distance, section 3, r 3^av-’ ft . . 0.4917 

Average radial distance, section 4, r 4^ a v> ^ ’ 0»3333 


Length of chordwise trailing section, j, ft 0.050 

Length of spanwise trailing section, b, ft 0.3333 

Thickness of trapezoid at coolant passage, Tg, ft 0.021 

Thickness of trapezoid at trailing edge, T^, ft 0.003 

Thickness of rectangle, T, ft . Q . . 0.010 

Temperature at blade root, T r , . F 330 

Thermal conductivity, kg 

Btu/(hr) (ft) (°F) 15-210 

Btu/(sec) (ft) (°F) . .... I 0.00417 -0.0583 

Cooling -air (in contact with rotor) temperature, T a , °F 0 

Heat-transfer coefficient between cooling air and rotor, h a 

Btu/(hr) (sq ft)(°F) 30 

Btu/(sec)(sq ft) (°F) 0.00874 

Average value in section 3 of ^± ) ^ !X x,'5 

Water, Btu/(hr) (ft) (°F) 30,360 

Btu/(sec) (ft) (°F) 8.433 

Ethylene glycol, Btu/(hr) (ft) (^F) 8305 

Btu/(sec) (ft) (°F) 2.307 

Kerosene, Btu/(hr) (ft) (°F) 6533 

Btu/(sec) (ft) (°F) 1.815 

Average value in section 4 of 

Water, Btu/(hr) (ft) (°F) 6107 

Btu/( sec) (ft) (°F) 1.696 

Ethylene glycol, Btu/(hr) (ft) (°F) 1656’ 

Btu/(sec) (ft) (°F) 0.460 

Kerosene, Btu/(hr) (ft) (°F) ....... .. 1318 

Btu/(sec) (ft) (°F) . 0.366 
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One -Dimensional Spanwise Temperature Distribution 

• The constants of integration were found by solving equations (5) 
to (10) for the specific examples stated; a blade thermal conductivity 
of 15 Btu/ ( hr) (ft) (°F) was used. Blade coolant passages extending to 
within l/l6 inch of the blade tip were considered in blades with spans 

ranging from 1-jjg- to 4^r inches, with water as coolant, and for effective 

gas temperature ranging from 2000° to 5000° F in order to study the 
effect of blade length on the temperature distribution. The temperature 
distributions were determined, by use of the calculated integration 
constants, from equations (1) to (4) and are shown in figure 6. The 
general trend shown is a nearly constant temperature about equal to the 
coolant temperature through the rotor, a sharp temperature increase 
through the rim and the base of the blade, a nearly constant tempera- 
ture (called the prevalent blade temperature) through the liquid-cooled 
part of the blade, and another sharp temperature increase to a value 
approaching the effective gas temperature in the uncooled part of the 
blade. From figure 6, it can readily be seen that the same prevalent 
blade temperature prevails through blades of various lengths cooled to 
within l/l6 inch of the blade tip. 

In order to study the effect of coolant-passage length on the tem- 
perature distribution, other calculations were made, with water as 
coolant and for an effective gas temperature of 2000° F, for a 4^ inch 

blade span but with blade coolant passages of various lengths. The 
calculations, again determined by use of equations (l) to (4), are 
plotted in figure 7; for short coolant passages high-temperature gra- 
dients exist throughout the blade whereas for long, coolant passages 
practically no temperature gradient exists throughout most of the 
blades. 

Calculations were also made for the same blades with 4-inch cool- 
ant passages and for an' effective gas temperature of 2000°. F with 
water, ethylene glycol, and kerosene considered as possible coolants. 
These results showing the effect of various coolants on temperature 
distribution are given in figure 8. For an effective gas temperature 
of 2000° F and a coolant flow of 7 pounds per minute per blade, the 
prevalent blade temperature in degrees Fahrenheit is about one fourth, 
two fifths, and one half of the effective gas temperature for water, 
ethylene glycol, and kerosene coolants, respectively. The consideration 
of -cooling air on the outside of the rotor accounts for the rotor tem- 
perature being less than the coolant temperature. 

Finally, for the spanwise case, figure 9 shows the variation of 
coo lant flow on prevalent blade temperature. The prevalent blade tem- 
perature decreases with increasing coolant flow; the rate of this 
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decrease diminishes as the coolant flow increases. It can be seen that 
the cooling effectiveness, defined as (Tg^ e -Tg)/(Tg e -Tj), changes « 

from 0.68 to 0.91 for water as the coolant flow is increased from 2 to 
16 pounds per minute per blade. The superiority of water is apparent. 


Three-Dimensional Temperature Distribution through 
Simplified Trailing Section 

A three-dimensional temperature distribution in a turbine -blade 
trailing section approximated by a rectangular parallelepiped (fig. 2) 
was determined by use of equation (12) (approximated by use of m = 1 
only, as previously mentioned) . The constants & n , <2 n , and were 

obtained from equations (16), (17), and (18), respectively, the con- 
stants an< l n from equations (14) and (15), respectively, 

and the constants an< l ^l,n from equations (13). The distribu- 

tion was found in two planes representing the maximum and minimum tem- 
peratures for the Z-axis; the first plane was located at the side of 
the cross section of the rectangular parallelepiped and the second 
plane was on the median plane through the rectangular parallelepiped.. 
The results are shown in figure 10(a) for an effective gas temperature 
of 2000° F and a thermal conductivity of 15 Biu/(hr) (ft) (°F) and give 
the temperature distribution at various distances from the coolant 
passage. The curve labeled "approximation of temperature" in fig- 
ure 10(a) is a one -dimensional chordwise distribution through the 
approximated rectangular trailing section. Similar results, for ther- 
mal conductivities of 120 and 210/Btu/(hr) (ft) (°F) are given in fig- 
ures 10(b) and 10(c),- respectively. Calculation of cos 
(see equation (12)) reveals the temperature variation in the two planes 
to be about 3.6 percent of Q for a blade thermal conductivity of 
15 Btu/(hr) (ft) (°F) , 0.5 percent for a thermal conductivity of 
120 Btu/(hr) (ft) (°F) , and 0.3 percent for a thermal conductivity of 
210 Btu/(hr) (ft) (°F) . For a thermal conductivity of 
15 Btu/(hr) (ft) (*%) , figure 10(a) shows a constant temperature for 
the last three-quarters of the blade; that is, conduction to the rim 
affects about the first quarter of the blade length. Figure 10(a) 
also shows that the level of the temperature in the principal portion 
of the blade rises rapidly as the distance from the coolant passage 
is increased. Near the. coolant passage, the distribution is in good • 
agreement with the one-dimensional spanwise distribution presented in 
figures. 6 and 8. As the thermal conductivity is increased, as shown 
in figures 10(b) and 10(c), the part ,of the blade affected by rim 
conduction increases; for a thermal conductivity of 

210 Btu/(hr) (ft) (°F) , about two thirds of the blade length shows this 
effect. 
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One -Dimens ional Chordwise Temperature Distribution 
through Simplified Shapes 

Because the three-dimensional distribution resulted in a con- 
stant spamri.se blade temperature in the region of the blade beyond the 
influence of rim cooling and because the critical blade point, as deter- 
mined from a temperature -stress relation, may likewise be beyond the 
region of rim-cooling influence, one -dimensional chordwise temperature 
distributions were obtained. In order to compare the three-dimensional 
distribution with a one -dimensional chordwise distribution, a rectan- 
gular trailing section was first considered.. Figure 11 shows this com- 
parison for an effective gas temperature of 2000° F and a blade thermal 
conductivity of 15 Btu/(hr) (ft) (°F) > the one -dimens ional distribution 
was determined by use of equation (19). From figure li it can be seen 
that the trailing-edge temperature increases from 1500° to 1850° F as 
the distance from the blade root is increased from 0.25 to 1.0 inch for 
the three-dimensional case, even for low-conductivity material. On the 
other hand, in the region of the blade beyond the influence of rim 
cooling (when the distance from the blade root is 1 in. or more), the 
one-dimensional and three-dimensional results approach identity; for 
such a region, a three-dimensional solution is unnecessary when low- 
conductivity materials are considered. 

In order to more nearly approximate a blade trailing section, a 
trapezoidal approximation was considered. For comparative purposes, 
the trapezoidal and rectangular sections were constructed to have equal 
lengths and areas. In general, however, the trapezoid is constructed 
so that the thickness at the coolant passages equals one -half the actual 
blade surface exposed to the coolant. This dimension, in turn, fixes 
the length of the trapezoid. 

The- temperature distribution through a trapezoidal section was 
obtained by use of equation (20) .’ Such a distribution, for an effec- 
tive gas temperature of 2000° F and a blade thermal conductivity of 
15 Btu/(hr) (ft) (°F) , is compared With the distribution through a rec- 
tangular section in figure 12. The temperature distribution for the 
trapezoidal section has a slightly steeper slope than that for the 
rectangular section at distances remote from the coolant passages and 
the temperatures at the trailing edge and. the coolant -pas sage wall are 
lower. Part of this lower temperature for the trapezoidal section is 
due to the additional thickness of the cooling surface. 

The effect of varying thermal conductivity from 15 to 
210 Btu/(hr) (ft) (°F) on the temperature distributions along the center ■ 
line of rectangular and trapezoidal sections is shown in figures 13(a) 
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and 13(b), respectively. In each case, as thermal conductivity is 
increased, the temperature distributions flatten and approach linearity, 
decreasing in the trailing section and increasing near the coolant 
passages. Figure 14 shows blade trailing-edge temperatures for various 
thermal conductivities obtained by use of the rectangular and trape- 
zoidal approximations. Cooling is substantially greater for the more 
representative trapezoidal section. For currently used high-temperature 
alloys with thermal conductivities in the neighborhood of 
15 Btu/(hr) (ft) (°F) , figure 14 shows a trailing-edge temperature dif- 
ference for the rectangular and trapezoidal sections of only about 
60° F for an effective gas temperature of 2000° F . 

Chordwise temperatures were also calculated, by use of equa- 
tions (22) to (27), for concentric circle annuli and parallel-plate 
approximations. These temperatures are not plotted in this report; 
they were only used as initial approximations for the two-dimensional 
numerical calculations made for an actual blade shape. 


Two-Dimensional Chordwise Temperature Distribution 

Two-dimensional temperature distributions were determined for an 
actual blade shape by application of the relaxation method. Separate 
calculations were made for the cases where an average gas-to-blade 
heat-transfer coefficient and where a typical variation in the gas- 
to-blade heat-transfer coefficient as shown in figure 15 were consid- 
ered. Initial trial solutions were determined by use of equation (20) 
for the trailing section, equations (22) to (24) for the leading sec- 
tion, and equations (26) and (27) for the thin-wall sections of the 
blade near the coolant passages. A comparison of the temperature dis- 
tributions for assumed variable and constant gas-to-blade heat-transfer 
coefficients, for an effective gas temperature of 2000° F, an average 
coolant temperature of 200° F, and a blade thermal conductivity of 
15 Btu/(hr) (ft) (°F) is shown in figure 16. The blade temperatures 
obtained are nearly the same for both cases except at the leading and 
trailing edges. Use of the average coefficient gives a conservative 
estimate of the trailing-edge temperature and a temperature that is 
somewhat too low near the leading edge. 

i 

A calculation has also been made for the average gas-to-blade 
coefficient and for a blade with a thermal conductivity of 
100 Btu/ (hr) (ft) (°F) . A comparison of this solution with, the similar 
one for a thermal conductivity of 15 Btu/(hr) (ft) (°F) is shown in fig- 
ure 17. The high-conductivity blade has about a 250° and a 600° F 
lower temperature than the low-conductivity blade at the leading and 
trai ling edges, respectively. Little temperature difference is obtained 
in the center of the blade, where extremely good cooling prevails. 
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The results of the various relaxation solutions are in good agree- 
ment with the one -dimensional chordwise approximations that were used 
to start the numerical solutions; , that is, good representative tem- 
peratures are obtainable by vise of one-dimensional chordwise approx- 
imations . 

In order to determine the effect of distance from the coolant 
passage on the trailing-edge temperature, trailing-edge temperatures 
were deter min ed for the liquid-cooled blade with five coolant passages 
shown in figure 18. Various length trailing sections were obtained by 
successively reducing the length of the trailing section shown in 
figure 18. The temperatures were determined by use of equation (20) 
and are shown for various thermal conductivities in figure 19. The 
trailing-edge temperature is reduced almost linearly as the length of 
the trailing section is decreased. The effect of thermal conductivity 
also decreases as the trailing -section length is decreased. 


NONDIMENSIONAL CHARTS 

The availability of several nondimens ional charts, to be subse- 
quently. discussed, eliminates the necessity for some numerical calcu- 
lation. 

The prevalent blade temperature is given by equation (2a) as 


m ho^oTg,e + hi,22i, 2^1 

T B,P ~ 


After division by hi^2H,2> 


1+X 


h o^o 
" " h i,2*i,2 

Subtraction of Tg^ e from both members of this equation leads to 

T g,e~ T B,p _ _1_ 

T g , e -T z 1+X 


B,p 

where 


hi + h. 0 
oo i ^ 2 i y u 




this equation may be written 


T 7 + AT 


e* . e 
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The plot of 
any given blade, 


T g,e- T B,p 

T g>e" T l 
X can be 


against X is shown in figure 20. For 

T -Tr 

evaluated and the value of can 

T g,e- T l 


be obtained. Finally, for the desired effective gas and coolant tem- 
peratures, a single simple algebraic operation results in the desired 
value of the prevalent blade temperature. 


Another nondimens ional chart, which gives the one -dimensional 
chordwise temperature distribution through a rectangular section, is 
also presented in figure 21. The temperature distribution through a 
rectangular section is given by equation (19), v 


T 


g> e 




i — (T P -T_) coshCpy' 
k B g> e * 


h i 

cp sinh j ' + cosh cpj * 
. B 


(19) 


For a given turbine and set of turbine operating conditions, values 
can be determined for all the quantities in this equation. A semilog 
plot with cpj ' as abscissa and hj_/cpkg as a parameter results in a 

nondimensional cooling ratio 


T 6,e- T B 


^ T g,e" T Z^ 


cosh cpy ' 


The addition of a second quadrant, with cosh cpy' as parameter, yields 

T g,e“ T B 

values of the temperature ratio — — at any point in the rectan- 


T g,e -T Z 

gular section. This 'chart is given in figure 21. The chart is used 
as follows: A vertical line is constructed through the calculated' 

value of the abscissa cp j * and is extended to the calculated value of 
the parameter hi/cpiq}. From this point, a horizontal line is drawn, 
extending into the second quadrant, and intersecting several lines 
representing various values of coshCpy 1 . Vertical lines from the 
intersection points to the abscissa in the second quadrant give values 


of the temperature ratio 


T g,e- T B 


gular section, 
obtainable. 


at various positions in the rectan- 
Values of the temperature Tg are then easily 


T g,e" T Z 
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In order to illustrate the use of figure 21, a recta ng ular section 
with dimensions equal to those previously tabulated will be considered. 
For a blade thermal conductivity of 15 Btu/(hr) (ft) (°F) and outside and 
inside heat-transfer coefficients of 222 and 2370 Btu/(hr)(sq ft)(°F), 
respectively, the parameters used in figure 21 are evaluated as follows: 



1 


2 ( 222 ) 

15(0.01)_ 


2 

(0.055) = 3 


h i _ 2570 _ „ 

cpk B “ (54.4) (15) ~ 


The point represented by these parameters is given by the symbol O 
in the first quadrant of figure 21) the symbol □ on figure 21 desig- 
nates the corresponding point for a thermal conductivity of 
210 Btu/(hr) (ft) (°F) . A horizontal line passing through the 
symbol O and intersecting the family of lines in the second quadrant 
of figure 21 gives as the abscissa of the second quadrant the values of 


T 

the temperature ratio 6 | >e 
gular section as follows : 



for various positions in the rectan- 


T g,e~ T B 

Tg,e-Tz 


0,075 

for 

o 

II 

9 

0.115 

for 

cpy« = 1 

0.285 

for 

Cpy' = 2 

0.75 

for 

cpy’ = 3 


For Tg^e = 2000° F and Tj = 200° F, it follows that the blade - 
temperature calculation results in the following values: 


T b = 1865 

for 

cpy' = 0 

( trailing -edge temperature ) 

= 1793 

for 

cpy' = 1 


= 1487 

for 

cpy' = 2 


= 650 

for 

to 

11 

>> 

8- 

(temperature at coolant passage) 


These results compare favorably with the calculated distribution shown 
on figure 12. 
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CONCLUDING REMARKS 

Analyses have been presented for obtaining spanwise temperature 
distributions near a coolant passage, chordwise distributions (for both 
approximated and actual shapes) in regions where conduction to the rim 
is inappreciable, and three-dimensional distributions for approximated 
blade shapes. Numerical examples based on specific blade configurations 
and heat-transfer coefficients available from unclassif ied sources have 
been presented. Although the analyses are exact, the numerical values 
used in the calculations may not necessarily be the same as those for 
cooled turbines. The numerical examples have been presented to indicate 
the range of applicability of the various analyses and to present the 
general nature of temperature distributions in liquid-cooled turbine 
blades. In the following paragraphs the more important results are 
reviewed to show their general guidance for design studies. 

The three-dimensional temperature distribution includes conduction 
to the rim, whereas the two-dimensional and one -dimensional chordwise. 
distributions do not. Consequently, the temperatures obtained from a 
three-dimensional investigation are less than those otherwise obtained. 
As a result, it is advisable to determine the simplified solutions 
first; if the temperatures that result are not excessively high, a 
three-dimensional investigation is unnecessary; whereas if the resulting 
temperatures are large, a three-dimensional investigation can be made. 

The two-dimensional solution just referred to is of, necessity a 
numerical solution (because of the varying boundary condition caused 
by the shape of the blade) . It has been shown by the calculations pre- 
viously presented that for uniform outside heat-transfer coefficients, 
simplified one -dimensional distributions were in excellent agreement 
with the two-dimensional relaxation solution. When variable outside 
heat-transfer coefficients were considered, the simplified solutions 
indicated optimistic results and the two-dimensional relaxation solution 
appears essential. Increasing blade thermal conductivity results in 
raising the cooling surface temperature and lowering the trailing-edge 
temperature and thus is a more uniform blade temperature. 

One-dimensional spanwise distributions proved to be valid near the 
coolant passages. The temperature of the cooled part of the blade 
(prevalent blade temperature) was independent of blade length. If the 
inside heat -transfer coefficient is increased, the prevalent blade tem- 
perature decreases; it was shown for a particular case that the prev- 
alent blade temperature for water as coolant was only about half that 
for kerosene as coolant. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, October 27, 1950. 
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APPENDIX A 
SYMBOLS 


The following symbols are used in this report: 


A 

b 

C. 

D 


area of metal , sq ft 
blade span, ft 

' ’ ^11 constants of integration 

diameter of circle used to approximate leading section 
of blade , ft 


d • 

H 1 

iH o 

h 



j 


ratio of distance between net points 0 and 5 and net 
spacing 

Hankel functions (special kinds of Bessel function) 

heat- transfer coefficient, Btu/(sec)(sq ft)(°F) or 
Btu/(hr)(sq ft) (°F) 

Bessel functions 


chordwise distance from blade trailing edge to 
coolant passage, ft 


K 

k 

l 


M,N,0,P 


Q 


1 



thermal conductivity, Btu/(sec) (ft) (°F) or 
Btu/ (hr) (ft) (°F) 

perimeter , ft 

points on figure 3 

heat flow, Btu/sec 
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r radial distance from hub of turbine, ft 
T temperature, °F 

x ' distance from blade tip to blade element, ft 
y distance from trailing edge to blade element, ft 
Z number of blades 

z distance from median plane of section to blade element, ft 

r H^i^Jod^) + iJlU^iHodSg) 


’ , T,(l-tan 

2k L^ + — - t - ^ -— j 


6 T g,e- T B 


0 

T 


hi ' 
o o 


h i,-2H,2 

distance between parallel plates, ft 
thickness of trailing section, ft 


cp 


2h^ 

V 


tan 


-1 / T 2" T 1 


2j 


Subscripts : 
a air 
av average 
B blade 


e 


effective 
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g gas 

i inside (inner with D) 

l liquid 



summation indices 


o outside (outer with D) 
p prevalent 

R rotor 

r blade root 

rim rim 

T blade tip 


1 1 


2 

3 


denote sections in spanwise investigation when used with x; 
1 and 2 denote trapezoidal thicknesses at trailing edge 
and coolant passage, respectively, when used with T 


4 J 

a 

denote end near blade tip and end near rotor hub of various 
sections • 

P . 


Superscript: 


prime linear dimension increased by T/2 
Functions : 


1 
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4h i f-1^ 11 " 1 

y> - frp _ m ^ V ■*•/ 

u g,e H } 2n-l 


/sin St-y 

2 ( T g,e' T B,r) ( g^y 


sin 2 Q n y 


,1 + 


2 V 
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APPENDIX B 

ONE-DIMENSIONAL SPANWISE TEMPERATURE DISTRIBUTION 

The heat balance for the uncooled section of the blade (fig. 1, 
section l) is given as follows: - • • . 

dT B i 

Heat entering by radial conduction = - k^A^ ^ 

Heat entering sides by convection = h Q Z 0 (Tg e -Tg j) dx-^ 


Heat. leaving by radial 


d ( d ^n t \ 

conduction = - ^ (^,1 * ^ j 


The heat-balance equation is 


dT. 


Vb,i + Vo\ e -Vi> “*1 - 

. V . jL( t + f!ki ax ) 

nn , 1 dx^ \ c , 1 + dx ] _ “l / 


or 


where 


d 2 T. 


dx. 


^-^ 2 T ' 2 ' 


B,1 


sd T 




stf - 



(Bl) 


A solution is 


T B,1 - T g,e - C 1 C0Sh ^< X 1 + C 2> 
where and Cg are integration constants . 


( 1 ) 
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For the cooled section of the blade (fig. 1, section 2) the heat 
balance is given as follows: 


dT- 


B,2 


Heat entering by radial conduction = - kgAg 2 ^ 

2 

Heat entering sides by convection = h 0 2 0 (T g ^ e -Tj^2) ^2 

d ( g \ 

Heat leaving by radial conduction = - dx 2 J 

Heat leaving by radial convection = h^ 2^ T B,2" T Z^ dx 2 


The heat-balance equation is 


dT, 


" dx* + h oV T g,e _T B,2^ ^2 

d ( dT B 2 \ 

^B^B,2 dxT 1 T B,2 + dXp ^2 J + h i,2 Z i,2^ T B,2" T Z^ **2 


or 


d 2 T- 


B , 2 


dx. 


^ 2 T. 


B,2 


= - 


(B2) 


where 


h o Z o + h i,2 l i,2 
^B^B,2 


and 


/ ^o^o T g,e + h i,2^i,2^2 \ 

l '^B A Bj 2 / 


/ 
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A solution is 


S&x. 9 - &fx 9 

TL _ = T C e c - C e ^ 

B,2 B,p 3 4 


where Cj and are integration constants and 

^ ~ + h.* oZj 


(2) 


T _ / 3r\ _ - o z o T g,e + h i,2*i,g T Z 

B,p-u; - Vo + h i, 2*1,2 


For the rim section (assumed constant area) of the. rotor (fig. 1, 
section 3) , curvature was neglected because the rim thickness was small 
in comparison with the rim radius. In addition, average values of 
h i,3^i,3 r 3 were used. The heat-balance equation reduced to 


d 2 T. 


rim 


dx„ 


- $ T. 


rim 


(B3) 


where 


and 


A solution is 


,2 _ ^ Jtr 5,av^ 1 a * ^i,5^i,3^av 
^rinAr im 


r 2 _ ^ rtr 3 , av^a^a + ^i,3^i,3^av^Z 
ArinArim 


<?x-z - Syi’z 

■'■rim = ^ + ^5 e + • ^6 e 

where Cg and Cg are integration constants and where 

2 


. ( 3 ) 
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For the rotor section (assumed constant strength) (fig. 1, 
section 4), average values are used for h^^J^ ^ and for r^_, and an 

approximate solution (as in section 3) is obtained. Because the rotor 
is liquid- cooled, such a solution is adequate. The resulting differ- 
ential equation for this section is 

2 

d Tp o o 

— - T R = - (B4) 


where 


and 


4nr 4 , av^a + ( h i,4 Z i,4)av 
^ t R A R 


2 _ 4jtr 4 , av^a-^a + ^i,4^i,4^ av-^Z 


r = 


k R A R 


A solution is 

Tr = / + c ? cosh (4) 

where Cy is an integration constant and where 

/■(/f 

The boundary condition 

dT 

- — = 0 for r = 0 


has already been applied. 
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APPENDIX C 

THREE-DIMENSIONAL TEMPERATURE DISTRIBUTION THROUGH 
SIMPLIFIED TRAILING SECTION 

A three-dimensional temperature distribution through a rectangu- 
lar parallelepiped, an approximation for the trailing section of a 
liquid-cooled turbine blade, is given as follows (fig. 2): 


ST- 


B 


Heat entering element from top = - kg dy dz — ^ 


ST- 


B 


Heat entering element from right end = - kg dx dz — ^ 


St. 


B 


Heat entering element from front = - kg dx dy 


_ a t b 


Heat leaving element at bottom = - dy dz ^ 


(STg S T b 


Heat leaving element at left end = - kg dx dz ^ ^ ^ dy 


'ST b S^Tg 

Heat leaving element at rear = - kg dx dy + 2 j 

Sz 


The heat-balance equation is 


St, 


dy dz " kB dz ~ V ta dy “sf 


St 


b 


St 


b 




°t b St b 

Sx 2 


- “b ay dz \"s + rr ^1 - *8 az * 


St b ^ % 


Sy" 


dy - 


'ST* S^t b 
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or 


.2 .2 2 

H 

* ft ™ V 




dy 


c>z‘ 


(ci) 


Simplification in applying boundary conditions results from the use of 
the effective gas temperature T„ „ instead of the blade temperature 

T„ as the reference temperature,, that is, by use of the substitution 

D 

9 = T g,e-% 


Equation (Cl) then becomes 

a 2 e 


a 2 0 


dx 


,2 


dy 


+ ^£= o 

- 2 Sz 2 


(11) 


The origin of the coordinates chosen is shown in figure 2; the 
plane z = 0 is the median plane of the rectangular parallelepiped 
and from considerations of symmetry no heat flows across this median 
plane. The boundary conditions to be applied are 


c>0 


c)0 _ , 

= 0 when 

x' = 0 

= 0 when 

w 

O 

II 

>> 

= 0 when 
oz 

z = 0 

T g,e- T B,r " hen 

x' = b 

■ h i< VV 8 > 

when y 


(C2) 

(C3) 

(C4) 

(C5) 

(C6) 


and 


k B M " V When z = T / 2 

A solution of equation (ll) satisfying the three boundary con- 
ditions' (C2) , (C3), and (C4) is 


(C7) 
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00 00 

n cos * n x ' cosh •^ r m,ny' cos > m z + 

m=l n=l 




fi) cosh @ x' cos 



^m,n 

m,n 


where 

X, 

X, xt, J/, 0 , 

& , and 

among 

these 

constants are - 





2 _ 2 

2 




*^m,n 

+■ y> 

’ (m 



2 2 

2 




#L w = y 

m,n m 

+ , 



m 


( 12 ) 


r (m^ n — 1,2,3, * • •) 


(13) 


The boundary condition expressed by equation (C5) leads to a 
determination of i £ . Substitution of equation (12) into equa- 
tion (C5) gives 


oo oo 


T -T_ = > > (X cos X b' cosh y' cos V z + 

g,e B,r • m,n n m,n J •'m 

m=l rt=l 

0 cosh & b' cos SI y' cos Jf z) 
m,n m,n n m s 


(C8) 


It can be seen that the first member in the double summation in equa- 
tion (12) vanishes when cos Jz? n x' = Oj. hence, if 

■ (” ' 

from equation (C8) it follows that 

00 oo 

^m,n cosl1 cos ^n^' cos - /< m z (C9) 

m=l n=l 

This is a Fourier development along y' and z. 
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The sums are equated term by term. Therefore, is determined as 
any solution of 


hi j ' 


tan ( a j ' ) = 


B 


, ' n *B*n 

The development of equation (C7) leads to a determination of jf . 


(17) 


m 


00 05 



m^X n?=l 


- k^HL „ cos x' cosh Jl y' (v sin + 

a run, n n m, n l m m 2 / 


? cosh w x' cos Q y' f V sin j/ $ 
m,n *^m, n *^n J l ' / m 2^ 


or 


oo co ^ 


= 2 _^ 2 _, ' h o(^m,n cos ^n x ’ cosh An,ny’ cos A £ + 
m=l n=l 

^m,n cosh ^nx' cos ^ny' cos *^m f ) 


(C13) 


tan V -J 
2 


k 

B^m 


V 

m"*" 


(18) 


The possibility of determining values of X m n and 0^,n t° 

satisfy equations (C9) and (Cll) has been established in textbooks on 
Fourier and other harmonic series (reference 7, pp. 118-121), and it 
is only necessary that the values determined define a convergent 
series. 


Values of are determined by integrating equation (Cll) 

m, n 

between the limits x' = 0 to x' = b' and z = 0 to z = t/ 2 and 
substituting the values previously determined for y n and in 

equations (16) and (18) . The integration is accomplished in two steps 
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using the functions cos 3^ x' dx' and cosj^ z dz 
where s and v are integers . For abbreviation, 

x = k B*^m,n sinh ^m,nJ' + h i cosh ^m,n<j 

Then 

f'x 1 =b ' 


rx 1 =t 

X C0S -V I 

m=l ^x'=0 


cos jg* x' dx 1 
n 


= h l( T g,e- T l) 


•'X ' =b 1 


/x'-O 


cos j^x' dx' 


because all the terms 


cos 3^ n x' cos i^ s x' dx' 


vanish if s ^ n. Integration of equation (Cl5) and 
limits leads to the result 


Z cos J'lz sin 2 i£. n x'"l^ 

^m,n 2i£ n ^n x ' + 2 J 0 

m=l ' ’ — 


= h i 




-ib' 


I J-.n 




sin 3? n x ' 


-10 


or 


00 




sin 2 3^_b' 

cos ^z l b * + 


m=l 


23? 

n 


sin y b' 
n 


£ 


n 


. as multipliers, 

(C14) 


(C15) 


substitution of 


(C16) 


- 2h i< T g,e- T I> 


(C17) 
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When equation (16) is used for’ <£ , there results 


00 n-1 

X **■«, n cos -V - ! V T S ,e'V Tifer ■ y n 

m=l 


When -equation (C18) is integrated in terms of z, using the 


multiplier cosjl^z dz, it follows that 




,ZnI 


0X, 


m ’ n J z =o 

because all the terms 


COS ^Kn Z dZ = 

m n 


f z= 0 


cos /z dz 


2 


cos ^z cos ^.z dz 


vanish if m ^ v. Upon integration, 


**n,nt „ v 2 


or 


^ ( T + 3111 2 2 \ 2 y n Sin/ 1 2 


m > n \ 2 ' 2^ 


X 


m 


Therefore, 


jr 


2^ n sin -^m 2 1 


^mi 


di/ sin 2^ IN ^ 

1 + — 

2 * y m 2 


(C18) 


(C19) 


(C20) 


(C21) 


(C22) 
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Similarly, from equation (C9) it is found that 


2 y- 
n 


sln ^ m 2 


cosh 0> b' 1 + 
m,n l 


sin- 2>J/ J 
m2 

2 ^ml 


(C23) 


where 


sin @ n j ' 

•^n = 2 ( T g,e _T B,r^ ~ sin 2 S n j' 

1 1 O * t 


2^ n j* 


(C24) 


Finally, substitution of equations (C14) and (C18) in equation (C22) and 
substitution of equation (C24) in equation (C23) lead to the following 
values of the coefficients and 


(-l) a_1 sin * yl m 2 8h i( T g,e" T l) 
r m,n = rt(2n-l) T „ T 

2 sin 2 ^m 2 

V + 2 ^m i 


^B^m,n s i r ^ 1 * / ^m,n^ ' + ^i cos k*An,nj~ 


sin^ n j' 4(T g>e -T B>r ) ^in 


'■ sin 2 g j ' \ 

1 + >m2 


sin 2 ^m 2 

2 >mi 


Tv cosh^L _b' 
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APPENDIX D 


ONE-DIMENSIONAL CHORDWISE TEMPERATURE DISTRIBUTION 
THROUGH SIMPLIFIED SHAPES 


Rectangular trailing section of blade removed from influence of 
rim cooling. - The heat "balance for an element of the rectangular 
section is given as follows (fig. 4(a)): 


Heat entering by conduction (right end) * kgT — 
Heat entering by convection (sides) = 2h Q 0 dy 


Heat leaving by conduction (left end) 




The heat-balance equation is 



or 


qp 2 0 


1 



A solution is 


d 2 0 

dy 2 


where 


o = C g cosh cp(y + C ) 


(Dl) 


(D2) 


where Cg and Cg are integration constants. The boundary con- 
ditions to be applied in the evaluation of the integration constants are 
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k B T § - h i T (V T l) - h i T ( T g,e- T Z- 0 ) when y =;J . 0»). 

and 

^ ^ = V when y ■ 0 (re) 

\ 

From equations (D2) and (D4) , it follows that 


k B < -'8 C f > sinh (cpCg) = cosh (CpCg) 

and therefore 

to*' 1 ® tanh-^ t • 

Cg * $ = y 2 (re) 

vheh T is not too large. Calculations revealed that the difference 
between Cg and T/2 varied from 0.4 to 0.5 percent for the values of 

cp considered and had no appreciable effect 0.005 foot or more inside 
the trailing edge. 

From equations (D2), (D3 ) , and (D5) it follows that 
kgTCgCpsinh (cpj 1 ) = V[ T g, e - T l-C 8 cosh 

and 'therefore 

h i 

fc( T g,e“ T l) (D6) 

C 8 “ ^ h— 

cp sinhf^j') +,t-— cosh (cpj ') 

*B 


Substitution of equations (D5) and (D6) in equation (D2) leads to the 
final equation 


4 (T g,e- T z) cosh ( cp y') 

0 „ j (19) 

cp si nh(cpj ') + ~ cosh (cpj ') 
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Trapezoidal trailing section of a blade removed from influence of 
rim cooling. - In reference 8, the temperature-distribution equation 
along the axis of a trapezoidal section is derived. The solution 
obtained is expressed by the relation 


e '= c io J o (i ^ + c n iH o (i0 (d?) 

where and C,, are integration constants, J and iH are 

XU XX O O 

Bessel functions, and £ is defined as 


where 



4K 


y* + 


T^(l - tan \|r) 

2 tan J 


(D8) 


and 



kg sin i|r 



tan 


V T 1 

2 j 


(D9) 


(DIO) 


The constants and in equation (D7) may be evaluated by 

application of the following boundary conditions: 

= 0 when y' = 0 (Dll) 

and 


- h i<Ve- T i- e ) r - J' 

or 


§£ = 0 when' £ ^ 


(D12) 


(D13) 


and 
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k 2K 2 

~d0 = h i( T g,e- T Z- 0 ) when C = ^2 

^2 a£ 


It follows from equation (D8) that when y 1 =0, 

[t, (1 - tan \jr)‘ 

«!-*! 

and when y' = j ' , 


2 tan \|r 


E 2 = at 


T^(l - tan \|r) 
J 2 tan i|r ‘ . 


1 

2 


Differentiation of equation (D8) gives 

d6 g d 9 
dT 2K 2 dy- 


hence, 


d0 _ . d0 

af = 0 vhen W ‘ 0 


From the properties of the Bessel functions, 

aj 0 (y) 


ay 


= - J x (y) 


and 


Mpfr) _ „ , , 

dy “ H i^ 


Differentiation of equation (D7) therefore gives 


(D14) 


(D15) 


(D16) 


(D17) 
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and it follows that 


when 


0 10 iJi(iJi) 

C ll- 


(D19) 


From equation (D12), with the use of equations (D7) and (D17) , it is 
found that 

2K^ 

k B XT [ _C 10 iJ l( i ^2) + c n H i( i ^ 2 )] 

u 

= ^[Tg^-Tl-CloJoU^) - Cni^dtg)] (D20) 

The values of the integration constants C^q anfl C . | 1 are now 

found by solving simultaneously equations (D19) and (D20) . Insertion 
of these values in equation (D7) gives the desired equation for the 
temperature distribution through a trapezoidal section 


0 = 


Vz. 

2 *% 


(T g,e- T l> K< 1 5 1 >j 0 < 1 » :+ 




- [H 1 d ?1 )iJ 1 d53)] + r 


2K % 


( 20 ) 


where 


For a wedge, T = 0. Then, from equation (D15) it is seen that 
= 0, and as a consequence, H-^.(i^) = ob and iJ^(i^) =0. 
Equation (20) then simplifies to 
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^■±^2 

— F- < T g,e- T l) 


( 21 ) 


Concentric circle annulus approximation for blade leading edge. - 
The derivations of the equations giving blade temperatures at the inner 
and outer edges of the annular region are given as follows (fig. 4(c)): 


o 

Heat entering from hot gas = 2n — h 0 ( T g,e“ T B,o) _(D2l) 

cLT 

Heat flowing over circle with diameter D = 2lt § _S (D22) 

d (i) 


D i 

Heat leaving to coolant = 2« -tj- h^Tg .j-Tj) 
A solution of equation (D22) is 



(D23) 


(D24) 


where C 12 is an integration constant. Application of the "boundary 
conditions 


and 


T B " T B,o 


^3 “ ^B,! 


when 


when 


D = D 0 
D = D. 

l 


the elimination of C 12 



to equation (D24) results in 
of Q, namely. 


and an evaluation 
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or 


„ 2rtk B( T B,o- T B,i) 

V D ° 
l0g e D ± 


(D25) 


From simultaneous solutions of equations (D21) (D23) , and (D25) it 
is found that 


Vo’Vi". kg 


T g,e-^ 


(25) 


1 + 


l0g e d 4 


Do (doV + D l h i) 


From a simultaneous solution of equation (D24) and equation (D24) with 
the inner "boundary condition applied and use of equation (D25 ) , it is 
found that 


log, 


e D. 


"*3 “ ^B,i + D ^ T B,o" T B,i^ 


(24) 


n o 

los e 


From equations (D23) and. (D25) and from equations (D2l) and (D25) 
there are obtained 


*B,i = T Z + 


2k B 


D q (T B,o _T B,i ) 
®i^i ^°®e D. 


( 22 ) 


and 


■V T g,e 


2k B ! 


Dh log„ — 
O O e D-: 


D 0 (T B,o" T B,l ) 


(23) 


Equations (22) to (25) express the blade temperatures in terms of known 
quantities. 
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Section of blade approximated by parallel plates. - The derivation 
of the equations expressing the hlade temperatures at the inner an^ 
outer edges of a blade section approximated by a region between parallel 
plates is given as follows (fig. 4(d)): 

Q = heat entering from the hot gas = h 0 (Tg^ e -Tg^ 0 ) 

Q = heat leaving to the coolant = ^(T-g^-Tj) 

K 

Q = heat flowing through the section = — (Tg 0 -Tg ^) 

where a is the distance between the plates. Equating these heat flows 
results in the following system of equations: 


ko 

~a~ (■%,o - -%,i) = ^o^g,e - -^B, o^ = ^i^B, i"^l^ 
A simultaneous solution of equations (D26) gives 

T Z + chh i ( h o T g,e +h i T l) 

rn _ u x 

1 + (h o +h i) 


and 


T g,e +> ah h ( h o T g,e +h i T z) 

T = °J: 

B,p k_. 

1 + -c=c- (h_+h^ ) 
al^h.^ v o i' 

From equations (26) and (27) it follows that 


T B,o“ T B,i ~ 


T g,e~ T l 


B 


1 + ^E7 <V*i> 


(D26) 


(26) 


(27) 


(D27) 
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Figure 1. 


- Arrangement of internal- 


cooling passages in turbine. 
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Figure 2. - Rectangular parallelepiped used to approximate 
trailing section of 4-lnch turbine blade for three- 
dimensional analysis. 



Figure 3. - Correction for heat received by trailing edge. 
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T g,e 


(a) Rectangle. 



(b) Trapezoid. 



Figure 4. - Simplified shapes used in cne-dimensional 

chordwise analysis . 
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Figure 6. - Effect of various blade lengths on temperature distribution for water 
cooled gas turbine for effective gas temperatures from 2000° to 5000° F. One- 
dimensional spanwlse analysis. Coolant passages extend to within 1/16 inch of 
blade tip; coolant flow, 7 pounds per minute per blade; thermal conductivity of 
blade, 15 Btu/(hr) ( ft) (°F) . 
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£ ‘ aanq-eaadiuaj, 


Figure 8. - Temperature distribution for gas turbine for various liquid coolants. 
One-dimensional spanwise analysis. Coolant-passage length, 4 inches; coolant 
passages extend to within 1/16 inch of blade tip; coolant flow, 7 pounds per 
minute per blade; thermal conductivity of blade, 15 Btu/(hr) (ft) (°F) . 
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(a) Thermal conductivity of blade, 

15 Btu/ ( hr ) ( f t ) ( °P ) . 

Figure 10. - Three-dimensional temperature distribution 
in trailing section of turbine blade. 
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(b) Thermal conductivity of blade, 

120 Btu/(hr ) ( ft ) (°F) . 

Figure 10. - Continued. Three-dimensional temperature 
distribution in trailing section of turbine blade. 
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(c) Thermal conductivity of blade, 

210 Btu/(hr) (ft) (°F) . 

Figure 10 » - Concluded. Three-dimensional temperature 
distribution in trailing section of turbine blade. 
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Figure 11. - Comparison of one-dimensional chordwise 
temperature distribution through rectangular trailing 
section with three-dimensional temperature distribution 
through rectangular parallelepiped. Thermal conduc- 
tivity of blade, 15 Btu/(hr) ( ft) (°F) . 



Temperature 
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Figure 12. - One-dimensional chordwise 
analysis of effect of shape on prevalent 
blade temperature. Thermal conductivity 
of blade, 15 Btu/(hr) (ft) (°F) . . 
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(a) Rectangular trailing section. 


Figure 13. - One-dimensional chordwise analysis of 
prevalent blade temperature for trailing section 
of blade. 
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(b) Trapezoidal trailing section. 


Figure 13. - Concluded. One -dimensional chordwise 
analysis of prevalent blade temperature for 
trailing section of blade. 
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Figure 14. - Trailing-edge temperatures obtained 
from one-dimensional chordwise analysis. Coolant 
temperature, 200° F. 
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gure 16. - Effect of variation of gas-to-blade heat-transfer 
coefficient on temperature distribution through cross section 
of water-cooled turbine blade. Thermal conductivity of blade. 
15 Btu/(hr ) ( f t) (°F) ; effective gas temperature, 2000° F; 
average water temperature, 200° F. 
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Figure 17. - Effect of variation of thermal conductivity on 
temperature distribution through water-cooled turbine blade. 
Gas-to-blade heat-transfer coefficient, 222 Btu/(hr) ( sq ft)(°F) 
effective gas temperature, 2000° F; average water temperature, 
200 ° F . 
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Length of trailing section of blade, in. 


Figure 19. - Trailing-edge temperature as function of 
length of trailing section of blade for water- 
cooled turbine blades with five coolant passages 
and various thermal conductivities. 
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Figure 20. - Nondimens ional chart for use in determination 
of prevalent blade temperature. 
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NACA-Langley - 4-18-51 - 1075 


Figure 21. - Nondimens ional chart for use in determination of temperature, distribution through rectangular trailing 

section of cooled turbine blade. 




One -dimensional spanwise temperature distributions One -dimensional spanwise temperature distributions 
gave satisfactory results near the coolant passages. gave satisfactory results near the coolant passages. 

One -dimensional chordwise distribution gave a good One -dimensional chordvn.se distribution gave a good 
first approximation to the actual solution in cases first approximation to the actual solution in cases 
where rim cooling was insignificant. -where rim cooling was insignificant. 
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Analytical methods are presented for computing 
:rature distributions in liquid-cooled turbine 





